Reactive glass and vegetation patterns 



(N 
O 
O 
(N 
-i— > 
O 

O 

in 



N.M. Shnerb 1 , P. Sara 2 ,H. Lavee 2 and S. Solomon 3 
1 Department of Physics, Judea and Samaria College, Ariel, Israel 44837 
2 Department of Geography, Bar-Ran University, Ramat-Gan, Israel 52900 and 
3 Racah institute of physics, The Hebrew University, Jerusalem Israel 91904 
(Dated: February 1, 2008) 

The formation of vegetation patterns in the arid and the semi-arid climatic zones is studied. 
Threshold for the biomass of the perennial flora is shown to be a relevant factor, leading to a frozen 
disordered patterns in the arid zone. In this "glassy" state, vegetation appears as a singular plant 
spots separated by irregular distances, and an indirect repulsive interaction among shrubs is induced 
by the competition for water. At higher precipitation rates, the diminish of hydrological losses in 
the presence of flora becomes important and yields spatial attraction and clustering of biomass. 
Turing-like patterns with characteristic length scale may emerge from the disordered structure due 
to this positive feedback instability. 

PACS numbers: 87.23.Cc, 89.75.Kd, 45.70.Qj 
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Vegetation patterns in the arid and the semi-arid cli- 
matic zones jl], |2] are an interesting example of sponta- 
neous symmetry breaking in complex systems. Compe- 
tition of shrubs for a limited supply of water is the rele- 
vant process that dictates the spatial organization. The 
struggle for water induces an indirect interaction among 
shrubs, as the flora extinct if its water supply is insuffi- 
cient. 

Competition for common resource has been considered 
for many years as one of the basic processes in population 
dynamics H [9. ft may be shown that, if two species 
compete for a common resource, the one that is able to 
survive at a lower resource level prevails and displaces 
the other species population. Stable coexistence of N 
species is possible if there exist N different resources and 
each of the species is a superior competitor for one of the 
supplies, that is, it has one biological niche. 

The situation becomes more complicated if the re- 
source admits spatial dynamics. Recent theoretical and 
experimental work reveals the dynamics of competing 
populations in water, where light, the limiting resource, 
is consumed gradually by the upper layers of aquatic phy- 
toplankton ||. This model may be extended to include 
spatial dynamics of the fauna, but it does not support 
time independent patterning . 

Vegetation patterns are an example of one species 
(shrubs) and one resource (water) system, where field 
studies revealed wide variety of stable, or almost stable, 
spontaneous segregation modes. Understanding the un- 
derlying mechanism for generation of such patterns and 
their observed resilience is considered as an important 
step toward a comprehension of the desertification pro- 
cess, where environmental effects like climate changing 
and grazing destroys the natural balance toward stable 
aridity. 

Technically, the water-biomass system has been con- 
sidered as a spatially extended nonlinear system, that, in 
some parameter range, may yield stripes, spots, labyrinth 
and other ordered arrangements attributed to a positive 
feedback mechanism, i.e., to the inhibition of water runoff 



and evaporation by the flora B. However, the typical 
perennial vegetation patterns in the semi-arid zone are 
disordered, as one can easily see in Fig. (I). The generic 
spatial organization of perennial flora varied along the 
precipitation gradient: from scattered "green spots" in 
the arid zone through clusters of shrubs in the semi-arid 
zone, to an almost full coverage of the soil by biomass in 
the humid/subhumid climate. Analysis of the transverse 
correlations in the three panels of Fig. (1) shows that 
the correlation length in the semi-arid zone is larger (by 
factor of 2-3) than in the other regions, and seems to in- 
dicate weak long-range oscillations in the Mediterranean 
site, perhaps a precursor of Turing instability. 

In this letter we present a general and simple model of 
the water-shrubs reaction that is able to yield all these 
features. Our model takes into account the intrinsic 
"noise", i.e., the amplification of initial fluctuations due 
to the minimal size needed for the survival of perennial 
flora. The resulting pattern is disordered but robust, thus 
it may be considered as the reactive equivalent of a glass. 
Although the concept of free energy, deep local minima 
and thermal equilibrium is absent in reactive systems, it 
still presents an example of a spontaneous breakdown of 
symmetry toward a disordered, long-living meta-stable 
state. 

To present the model, let us begin from its zero dimen- 
sional ("flower pot") deterministic and continuous dy- 
namics. With water supplied to the system at some rate 
R and continuum vegetation growth, the time evolution 
of the water-shrub system is described by the following, 
nondimensionalized rate equations: 



dB 

~dt 

dw 



wB - fiB 

R — w — XwB 



(1) 



Where w stands for the available water density, B is 
the density of shrubs biomass, the term wB represent 
shrubs growth as they consume water while —XwB is the 
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FIG. 1: Results of direct field measurements at three different locations along the precipitation gradient. The distribution of 
perennial shrubs (annual flora not included) is presented for an area of 100 square meters at each site. Each black spot presents 
a shrub and the size of a spot is proportional to the area size of the canopy. Shrubs distribution on hillslopes has been taken 
at three sites representing mildly-arid, semi-arid and subhumid climate conditions in Israel. The mildly-arid site (left panel, 
with mean annual rainfall 260mm, hillslope gradient 11°) is located at Mishor Adumim, 10km east of Jerusalem. The semi-arid 
site (middle, with mean annual rainfall 330mm and gradient of 16°), is located at Ma'ale Adumim, 8km east of Jerusalem. 
The Mediterranean site (right, mean annual rainfall 620mm, hillslope gradient 13°), is located at Giv'at Ye'arim, 11km west 
of Jerusalem. All the three sites have hard calcareous bedrock and southeast exposure (azimuth 140° — 150°). 



corresponding consumption of water by shrubs, (i is the 
"death rate" of the vegetation in the absence of water 
and the term — w represents water losses by percolation 
and evaporation. 

The set of differential equations (Q) admits two non- 
negative fixed points. The trivial one at Bq = 0, wq = R, 
becomes unstable to small perturbations as R = fj,, while 
above this value a "coexistence" fixed point at 

w 1 = f i Bi = $f£, (2) 

becomes a stable node (if R 2 < 4(R — ^<)/i 2 ) or a stable 
focus (if R 2 > 4(R- /i)/i 2 ). 

Adding lateral water flow to the above model leads to 
a reaction-diffusion equations of the form, 



— = wB-^B (3) 
duu 

— = DV 2 w + v-Vw + R-w- XwB, 
at 

where v pointed down the hillslope. Simple linear anal- 
ysis implies that in the absence of cross-diffusion effects 
(like those considered recently by Q), no Turing-like in- 
stability exists in that system; the steady state is a uni- 
form covering of all the plane by the same amount of flora 
which corresponds to the stable fixed point and fluctua- 
tions of wavenumber k decay like e~ k *. 

In the desert area considered here there are two sea- 
sons, dry summer and humid winter. Eqs. (§,§) presents 
the winter, with "smeared" rain events. While annual 
flora wilt in the summer, perennial shrubs have to sur- 
vive, so they must reach some threshold size before the 
dry season. If the winter is not long enough to allow for 
a full development of the plant to its stable fixed point, 



the survival of a shrub depends on its size at the end of 
the rainy season, which, in turn, depends on small fluc- 
tuations in its initial size and the consumption of water 
by its neighbors. 

In the next winter, the existence of a shrub causes a 
depletion of the available soil moisture in its immediate 
neighborhood (roughly speaking, in an area of typical lin- 
ear size and the chance for another shrub that 
pops out in the depletion region to reach the threshold is 
lowered. The whole area is then segregated into a mosaic 
of water accepting and water contributing patches. This 
is a non- Turing mechanism that has nothing with the 
effect of shrubs on the overland flow. While Turing in- 
stability is characterized by some typical wavelength that 
sets the linear size of vegetation and bare soil patches to 
be equal, our model allows for clusters of arbitrary size, 
as indicated in Fig. (1). 

The optimal segregation of the hillslope, that gives 
maximal biomass per unit area, is an ordered array of 
shrubs, each located in the lower end (or, on a flat plane, 
in the middle) of its contributing area. A regular or dis- 
torted "lattice" of flora is then formed, similar to the 
structure of atoms in a two dimensional crystal. How- 
ever, this optimal scenario is rarely accomplished in na- 
ture, due to the stochastic character of the growth pro- 
cess itself. For simplicity, let us assume that the seed 
bank in the soil ensures the development of a perennial 
shrub if some water exist at the site. As the first shrub 
pops out in an empty region, the soil moisture in its sur- 
roundings (primarily downhill) is depleted, and the next 
shrub will not grow in this "shadowed" area. Neverthe- 
less, as the next shrub also occurs at random, its position 
will be uncorrelated with the first, except that it can not 
pop out at the shadowed region of the first one, and so 
on. The process continues until all the slope is shaded 
( "jammed" ). This stochastic growth process yields a ran- 



FIG. 2: Numerical results of forward Euler integration of the 
reaction-diffusion equations (^) on 100X100 sites grid with 
periodic boundary conditions. A grey spot is plotted around 
the location of each shrub, with the size of the spot propor- 
tional to its biomass. The dotted lines present soil moisture 
contours around the shrubs. The simulation parameters are 
v = (no slope) and D = 10, /J, = 0.2, R = 0.5, A = 1.2 (this 
implies Bi = 1.25 and w\ = 0.2). Initial conditions are no 
water and a seed of biomass taken from a square distribution 
between [0,0.01]. Effect of the summer is modelled at the end 
of each "winter" (21 time cycles) by setting all water to zero 
while the flora at a site is dropped to the seed level if the 
biomass is smaller than a threshold, with B t h — 1.2. The av- 
erage (per site) values of water (about 0.32) and flora (0.57) 
reflect the inefficient use of water attributed to the glass-like 
structure. 



dom covering of the slope by shrubs, with typical distance 
between nearest neighbors, but no long-range structure. 
The random arrangement is, however, extremely robust; 
although the slope is covered inefficiently by the shrubs, 
there is not enough source area for the next shrub. Es- 
sential changes, such as the death of a plant and a for- 
mation of another one, are discontinuous, and though are 
very rare, unless some intervene come from the outside, 
perhaps in the form of grazing or climate changes. 

In Figure @ , typical results of a numerical simulation 
of a threshold-noise model on a "flat" (no slope) land are 
presented. In the simulation, the system freezes rapidly 
(in 1-2 rainy seasons) to a robust configuration (slower 
dynamics will increase the efficiency of biomass growth) , 
and then persists, with negligible fluctuations, up to 500 
winters. The average amount of flora is much smaller 
than Bi , and the average amount of water is larger than 
w±, i.e., there is an inefficient use of the water due to 
the stochastic arrangement of the shrubs. To guide the 
eye, soil moisture contours are also plotted (dotted lines, 
with the water level inside smaller than 0.46) and reveal 
the depletion zone around each shrub . No empty site 
in the region maintains enough water to allow for new 
shrub to develop, and the whole region is "shaded" by 
the existing flora. Figure (||) present the results of a sim- 
ulation with the same growth and diffusion parameters, 



FIG. 3: Same as Fig. (|3j) but with a downhill drift of water 
where the asymmetry parameter is v = 0.5. Note the change 
in the shape of the water contours. Average biomass and 
water levels are almost the same as in Figure (0). 



but with nonzero downhill slope, and the effect on the 
moisture depletion zones is evident. 

Various aspects of this competition scenario are simi- 
lar to the adsorption of large particles at an interface || . 
In the model of random sequential adsorption, "hard" 
particles are added sequentially to a D dimensional vol- 
ume at random positions with the condition that no trial 
particle can overlap previously inserted ones. The ad- 
dition process is than repeated until the system reaches 
its "jamming limit", at which the density saturates. The 
adsorbed particles density at the jamming limit is lower 
than the close-packed form, and the configuration of ad- 
sorbed particles is " freezed" at some disordered pattern. 
A shrub above threshold, with its excluded volume of 
depleted moisture is similar to an adsorbed "disc". The 
shrubs-water model, although non-local (excess water 
transferred downhill) and reversible (new shrub may re- 
move an existing plant by depleting its water resources), 
yields a similar jamming disordered and inefficient cov- 
ering of the slope. 

The above considerations about glassy structure eluci- 
dates the existence and the robustness of the vegetation 
patterns in the arid zone, but fails to explain the aggrega- 
tion of shrubs (larger correlation length) in the semi-arid 
zone, or the Turing patterns HQ]. To explain patchiness 
one should consider the positive-feedback mechanism, 
i.e., the inhibition of water dynamics induced by the 
shrubs themselves. In the absence of shrubs (and other 
meso topography factors) , water flow downhill, with some 
typical lateral displacement per unit length. In the pres- 
ence of shrubs the flow in their vicinity is slower than 
that in bare soil, as a result of higher infiltration rates. 
In addition the microclimate under shrub is character- 
ized by less direct radiation and smaller evaporation rate 
||. This means that close to the shrub there are favor- 
able soil water conditions and more flora may grow. The 
"repulsive" interaction among shrubs due to the struggle 
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FIG. 4: Effect of the positive feedback mechanism. While at low precipitation (R — 0.45, left panel) the repulsive interaction 
wins and no clustering occurs, at higher humidity (R=0.47, middle panel) clustering is more pronounced, and at R = 0.5 (right) 
vegetation stripes may be recognized. All other parameters are the same as in Figure (0). 



for water is then balanced by an "attraction". Accord- 
ingly, the size of a typical cluster changes along the pre- 
cipitation gradient, from a single shrub at the arid limit 
to large clusters in the semi-arid regime. The response 
of the system to external parameters (climate change) 
and grazing seems to depend on its phase. In particular, 
it seems that hysteresis loops (desertification transition) 
like those described in fij] are a characteristic of the clus- 
tered phase. A detailed discussion of these issues will be 
presented elsewhere. 

In Fig. (|j), the results of the same simulation program 
with positive feedback are presented. The only new in- 
gredient added to the simulation is a suppression of the 
asymmetry in the downhill water flow, with the asym- 
metry term v multiplied by exp[—5B(r) / B t h] (this is a 
cross-convection effect). Diffusion and evaporation re- 
main the same as in (0). As the downhill flow becomes 
smaller, water tend to accumulate at shrub's neighbor- 
hood. The transition from the arid (left) zone, with no 
clustering and glassy structure, to more clustered pat- 



terns and even some linear order in the semi-arid (right) 
is evident. 

In conclusion, the generic spatial patterns due to the 
struggle for water are disordered frozen patterns, and the 
threshold-noise process dictates the vegetation spatial or- 
ganization in the arid zone. The instability that yields 
this glassy structure is not Turing-like, thus the inho- 
mogeneity is not characterized by typical length scale. 
Shrub clusters and Turing patterns emerge as further in- 
stability of this glassy structure if the precipitation is 
large enough, where the positive feedback dominates. It 
then leads first to clustering of shrubs and then to a 
global order in the form of "tiger bush" patterns. 



Acknowledgments 

We wish to acknowledge E. Meron, F.J. Weissing and 
J. Huisman for most helpful comments and discussions. 



[1] D.J. Tongway anf J. A. Ludwig, Australian Journal of 

Ecology, 15, 23 (1990). 
[2] C. Valentin, J.M. Herbes and J. Poesen, Catena 37, 1 

(1999). 

[3] G. F. Gause, The Struggle for Existence, (Williams and 
Wilkins, Baltimore, 1934). 

[4] D. Tilman, Resource Competition and Community Struc- 
ture, (Princeton University Press, Princeton, 1982); J. 
D. Murray, Mathematical Biology (Springer- Verlag, New 
York, 1993). 

[5] J. Huisman and F. J. Weissing, Ecology 75 507 (1994); F. 
J. Weissing and J. Huisman, Jour, of Theoretical Biology, 
168, 323 (1994); J. Huisman, R.R. Jonker C. Zonneveld 
and F.J. Weissing, Ecology 80, 211 (1999). See also S. 
V. Petrovskii and H. Malchow, Theo. Pop. Bio. 59, 157 
(2001) for a prey-predator model of zooplankton that sup- 
ports oscillating spatial patterns. 

[6] J.B. Wilson and A.D.Q. Agnew, Adv. Ecol. Res. 23, 263 



(1992). R. Lefever and O. Lejeune, Bull. Math. Biol. 59, 
263 (1997). H. Lavee, A.C. Imerson and P. Sara, Land 
Degradation and development 9, 407 (1998). 

[7] J. von Hardenberg, E. Meron, S. Shachak and Y. Zarmi, 
Phys. Rev. Lett. 87 198101 (2001). 

[8] Actual field experiments usually finds the soil moisture 
under a shrub larger than the moisture few meters away 
from it. The reason for that has to do with the positive 
feedback mechanisms that reduce water losses in this re- 
gion. However, up to some limit, this water resource serves 
the shrub during the summer, and is not available for the 
use of another planet. One may regard this fraction of soil 
moisture as part of the biomass, while the contours in Figs. 
(QJ^Q) represent available water density. 

[9] See, e.g., J. Talbot, G. Tarjus, PR. Van-Tassel and P. 
Viot, Physica A 165, 287 (2000) and references therein. 



